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: Abstract 



' Nuclear parton distributions are studied in a parton model with rescaling 

ON ' 

' and recombination mechanisms. Parton x distributions are first calculated 
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00 

ON ' be compared with various F2 data. The experimental shadowing is explained 



at Q 2 =l GeV 2 by the model, and they are evolved to larger Q 2 in order to 



43 

Q_i' by the parton-recombination effects. Then, the modification of sea-quark and 

gluon distributions is investigated. We show x and Q 2 dependent results of 



the ratios Sa/Sd and Ga/Gd- They indicate significant shadowing at small 



^ ' x and large ratios in the medium-x region. Our theoretical results should be 



tested by future experiments, for example, at RHIC (Relativistic Heavy Ion 

Collider). Furthermore, our studies are important for finding a signature of 

the quark-gluon plasma in heavy- ion reactions. We also indicate that nuclear 

effects in the deuteron should be taken into account for finding the accurate 

F2 structure function of the neutron. 
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I. INTRODUCTION 



Although unpolarized parton distributions in the nucleon are now known in detail from 
very small x to relatively large x, the distributions in nuclei are not investigated to such a 
great extent [[Ij]. It was usually assumed in the 1970's that the nuclear structure functions 
F 2 were simply given by the one for the nucleon multiplied by the mass number for isoscalar 
nuclei: F 2 = AF 2 N . Nuclear binding energies are very small compared with energy scales 
in deep inelastic scattering, so that they were expected to be insignificant in discussing 
the structure functions. It is, therefore, rather startling to find significant modification 
of the nuclear F 2 at medium x in the European Muon Collaboration (EMC) data |2j in 
1983. The nuclear modification is now named the ("old") EMC effect after the experimental 
collaboration name. 

In the middle of 1980's, theoretical models were proposed for explaining the nuclear 
modification in the medium-x region. The first idea after the EMC finding is a Q 2 rescaling 
model || in 1983. The original motivation for the rescaling model was to consider that 
the nucleon size is modified in the nuclear medium. Because the structure function F 2 A 
reflects lightcone momentum distributions of nuclear partons, the size change results in the 
modification of F 2 . Although such modification may exist, the rescaling model is no more 
interpreted in the original form. It is better to be considered as an effective model [|J, 
which could be consistent with the following binding model. A conservative idea in terms of 
nuclear binding was, on the other hand, proposed in 1985 [f| . The nuclear F 2 is given by the 
nucleon's F 2 convoluted with a spectral function, which indicates momentum distributions of 
the nucleons in the nucleus. Due to the binding energy e, the x distribution of F 2 is shifted a 
few percent (~ e/m N ). Although the magnitude itself is not large, the shift becomes 10—20 
% effects in the F 2 at medium x. Because the model is successful in explaining a significant 
portion of the experimental data, the binding model is now accepted as the most probable 
interpretation of the EMC effect at medium x. The model can also explain the large-x 
region, where the ratio F 2 A /F 2 D increases as x — > 1, in terms of the nucleon Fermi motion. 
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We discuss the compatibility of the rescaling model with the binding model in section ||. 

The small- £ region was studied particularly from the late 1980's. Accurate ratios F^/F 2 D 
were obtained experimentally at small x for several nuclei by the New Muon Collaboration 
(NMC) || and the Fermilab-E665 collaboration jTJ. It became possible to test theoretical 
calculations in comparison with the data. The nuclear modification at small x is called 
shadowing, which means that internal constituents are shadowed due to the existence of 
nuclear surface constituents. A laboratory frame picture for describing the shadowing is a 
vector- meson-dominance (VMD) type model [§]. A virtual photon transforms into vector 
meson states, which then interact with a target nucleus. The propagation length of the 
hadronic fluctuation exceeds the average nucleon separation (2 fm) in nuclei at x < 0.1, so 
that the shadowing occurs because of multiple scattering of the hadron inside the nucleus. On 
the other hand, the shadowing phenomena could be explained by a parton-recombination 
model PHT2[, which is a description in an infinite momentum frame. The recombination 



occurs because the localization size of a parton with momentum fraction x exceeds 2 fm at 
x < 0.1. 

We showed in the parton model with the rescaling and recombination mechanisms that 
the nuclear F 2 could be explained from very small x to large x JTT|jr3| . Although the nu- 



clear F 2 structure functions are established, the sea-quark and gluon distributions are not 



well tested [14]]. Of course, because the F 2 at small x is essentially given by the sea-quark 
distribution, the sea modification should be equivalent to the F 2 modification in the leading 
order (LO) of a s . However, we would rather have independent confirmation, for example, 



by the Drell-Yan experiment. Although there exist the Fermilab-E772 Drell-Yan data |T5| 



the small-x shadowing part and the medium-x region were not probed. Furthermore, there 
is little experimental information on the nuclear gluon distributions pl6||L7| . An item of 
good news is that the Relativistic Heavy Ion Collider (RHIC) will be completed in the near 
future, so that the shadowing in the sea-quark and gluon distributions could be observed ex- 
perimentally by the Drell-Yan and direct-photon processes. Considering these experimental 
possibilities, we should try to predict the shadowing behavior theoretically in order to test 



our theoretical understanding on high-energy nuclear structure. 



The major purpose of this paper is extend the parton model in Ref. |TT[ to the studies of 
the sea-quark and gluon distributions in nuclei. In addition, the parton distributions in the 
nucleon are updated to the new ones. The small- x distributions were particularly modified 
due to the HERA F 2 data after the publication of Ref. [IT]. Another purpose is to study 
the parton distributions in the deuteron. They were assumed to be the same as those of the 
nucleon in the previous publications |ll~l , |T6| . In section |T|, our parton model is introduced 
with explanation for the Q 2 rescaling and parton-recombination mechanisms. In section 



llj, the F2 structure functions are calculated for various nuclei in order to be compared 
with measured experimental data. In section [IV], the model is extended to the sea-quark 
and gluon distributions. Their x and Q 2 dependent results are shown. As a byproduct of 
our investigation, we can discuss the relation between the deuteron and neutron structure 
functions. In section [V], our results are compared with the experimental ratios F^/F^- The 
summary is given in section [VT|. 



II. NUCLEAR MODEL 

People studied possible physics mechanisms for interpreting the medium-x modification 
of the F 2 structure functions just after the EMC finding. A few years later, the small- 
x mechanisms were investigated independently from the medium-x physics. It is rather 
unfortunate that enough efforts are not made to unify these pictures to produce F2 from 
small x to large x. For example, the parton distributions in the wide x range are necessary for 
studying Q 2 dependence of the structure functions because they should be supplied as input 
distributions. For this reason, we set up a parton model which could describe the parton 
distributions in the whole x region in a dynamically consistent way |Tl[|. We calculated the 
modification of the nuclear parton distributions at certain Q 2 (= Qq) by the rescaling and 
recombination mechanisms. Then, they are evolved to larger Q 2 where experimental data 
are taken. Because the same model is used in our studies, we explain only the essential part 
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in the following. The interested reader may look at the original papers for finding 

the detailed formalism. 



A. Q 2 rescaling 

As explained in the introduction, the Q 2 rescaling model was originally proposed as a 
first explicit quark signature in nuclear physics in terms of nucleon-size modification in the 
nuclear medium |J. As a result, the nuclear F 2 structure function is related to the nucleon's 
F 2 by a simple Q 2 rescaling, 

F 2 A (x,Q 2 ) = F 2 N (x^ A Q 2 ) with C A = (X 2 A /X 2 N r^ )MQ2) , (2.1) 

where £ A is called rescaling parameter. The factors A^ and A # are confinement radii for the 
nucleus A and the nucleon, respectively. 

The proponents of this model modified the interpretation ^fl. They rather consider the 
rescaling model as an effective model, which could include the binding and/or nucleon mod- 
ification effects. In order to justify this idea, they used factorization-scale independence. 
It means that physics observables should not depend on the scale which is introduced ar- 
tificially by hand. This is generally used in structure functions, which are separated into 
nonperturbative quantities (parton distributions) and perturbative ones (coefficient func- 
tions) by the separation scale /i. The scale fi is introduced artificially, so that the structure 
functions should be independent of the scale. It leads to so called renormalization group 
equation. In our situation, the factorization means that a nuclear structure function is 
separated into a nucleon momentum distribution and the corresponding structure function 
in the nucleon. The factorization scale independence is not, to be exact, proved in this 
case. However, we consider it rather as a working hypothesis: the observables, nuclear 
structure functions, are independent of the artificial separation into nuclear physics and 
subnucleon physics. Then, the nuclear binding model could be also expressed in term of the 
Q 2 rescaling terminology. For example, the binding energy and Fermi momentum in the 
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Fermi gas model are related to the rescaling parameter by || e/rriN = —■y^'^A/i^Po) and 
k 2 F /m 2 N = 5(2 T f s - 7f s )« A /(2A)) + 0(e 2 /m 2 N ), where 1 - k a = a s (UQ 2 )/a s (Q 2 )- There- 
fore, the rescaling model could be considered as an effective model which could include the 
binding-type nuclear effects and possibly also subnucleon-type contributions. In this sense, 
the rescaling parameter £4 could be determined by fitting some experimental data, rather 
than calculating it in a particular nuclear model such as the one in the rescaling paper of 
1985 |§. In any case, the original fits are no longer valid if Fermi-motion effects are included 
lilj . As explained in section [II, the parameter for the calcium nucleus £ca is determined 
by the / F® data at medium x. The other nuclear values are calculated by using the A 
dependence relation of the rescaling model. 

B. Parton recombination 

We use an infinite momentum picture, the parton- recombination model |9HT2|], for de- 
scribing the small- 2 shadowing physics, although there are other descriptions, for example, 



VMD-type models M and Pomeron models [19]. These different models explain the nuclear 



F 2 structure functions in some ways or other, so the appropriate description may be de- 
termined by studying other observables such as nuclear gluon distributions or possibly the 



valence-quark distributions at small x pO 



The parton recombination, which means the interaction of partons in different nucleons, 
is used as a model for interpreting the nuclear shadowing f9Hl2| . This mechanism could occur 



by the following reason. In an infinite momentum frame, the average longitudinal nucleon 
separation in a Lorentz contracted nucleus is L ~ (2.2 fm)M A / Pa = (2.2 fm)m N /p N , and 
the longitudinal localization size of a parton with momentum xp N is AL » l/(xp N ). At 
small x (x < 0.1), the localization size becomes smaller than the average nucleon separation. 
It means that partons from different nucleons could interact significantly. The interaction is 
called parton recombination or parton fusion. According to this mechanism, the modification 
of a parton distribution Ps(x 3 ) is given by 
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Ap 3 (x 3 ) =K J dx 1 dx 2 pi(x 1 ) p 2 (x 2 )T PlP2 ^ m (x 1 , x 2 , x 3 ) 5(x 3 - x x - x 2 ) , (2.2) 

for the recombination process p±p 2 — > p 3 . Here, the factor K is defined by K = 
^A 1 ' z a s /{2R 2 Q 2 ) with the expression for the nuclear radius R = -Ro^ 1 ^ 3 - The function 
T PlP2 ^ P3 (xi, x 2 , x 3 ) is called parton fusion function. Because the recombination is opposite 
to the splitting process, the T PlP2 ^ P3 is related to the corresponding splitting function P P1 <- P3 
in the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) Q 2 evolution equations: 

^PlP2~ J >P3 5 *^2) ^3) ~ Ppi^P3 (• r l/- r 3) CpiP2^P3 > (2-3) 

X 3 

where C PlP2 ^ PA is the color-factor ratio, for example C q c-, q = X! i^tiT^ti/ (tfkYtfk- The 

(l,a),k (k),l,a 

averages are taken over the initial color indices, for example, (k) indicates the average over 
the index k. There is also a momentum-cutoff factor w(x) = exp{—m 2 N ZQX 2 /2) for the leak- 
out partons. It indicates that large-momentum partons are confined in the small region so 
that they do not leak out from the nucleon. It is known that z$=2 fm is an appropriate value 
il|,|TT}1 , so that it is used throughout our analysis. Because the recombination formalism 



itself is discussed in detail in Refs. jlO|,[ll| , we do not repeat it here. The explicit equations 



are given in the Appendix part of Ref . JIT ] . 



III. STRUCTURE FUNCTIONS 

In order to calculate nuclear parton distributions, it is necessary to employ a set of parton 
distributions in the nucleon. The optimum parton distributions are obtained by fitting many 
experimental data from various processes. There are three major groups in the studies of 
the unpolarized parametrization. They are CTEQ (Coordinated Theoretical/Experimental 
Project on QCD Phenomenology and Tests of the Standard Model), GRV (Gliick, Reya, 
and Vogt), and MRS (Martin, Roberts, and Stirling). Among these parametrizations, we 
decided to take the MRS-R2 distributions |22[ without any serious reason. Of course, the 



details of our numerical results depend on the input distributions; however, the essential 
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results should not be altered. For simplicity, the MRS-R2 is modified to have the flavor 
symmetric sea, u = d = s = Sea/Q, in the following analyses. 

We have to determine the point Qq , at which the initial nuclear distributions are calcu- 
lated. It is considered as a parameter in our model, and it is determined in the following 
way. In order to explain the large shadowing, for example in F^/F®, it is necessary to take 
rather small Qq because the recombination effect has a higher-twist nature (~ 1/Qq). On 
the other hand, perturbative Quantum Chromodynamics (QCD) is assumed in evolving the 
calculated distributions at Ql to those at experimental Q 2 points. Therefore, Qq should be 
in the perturbative region. With these two reasons, we end up taking Qq ~1 GeV 2 . Because 
the MRS-R2 distributions are set up at Q 2 =l GeV 2 , it is simply taken as the MRS-R2 point, 
namely Qq=1 GeV 2 , in the following analyses. 

Using the MRS-R2 distributions, we calculate first the rescaling effects for the valence 
quark distributions at Qq =1 GeV 2 . Because this procedure violates the momentum con- 
servation, the gluon distribution is multiplied by a constant amount so as to satisfy the 
conservation. Then, the obtained distributions are used as input distributions for calculat- 
ing the recombination effects. The experimental Q 2 values are different depending on the 
x region: they are typically a few GeV 2 at small x and a few dozen GeV 2 at large x. In 
order to compare with the experimental data, we evolve the obtained nuclear distributions 
at Q 2 =l GeV 2 to those at Q 2 =5 GeV 2 by the parton- recombination (PR) Q 2 evolution 



equations |2l||24]j . In particular, we use the equations proposed by Mueller and Qiu p4j an d 
their numerical solution in Ref. |25]. The QCD scale parameter of the MRS-R2 is A=0.344 
GeV. 

In calculating the rescaling, it is necessary to find the rescaling parameter £4. As we men- 
tioned in the previous section, £c a is determined for explaining the experimental medium-x 
ratios, F 2 Ca /F 2 D , reasonably well. It is found that Cca— 1-74 is the appropriate value at Q 2 =l 
GeV 2 . The rescaling parameters of other nuclei are determined by the relationship of the 
original rescaling model. Namely, the scale parameter is determined by \ia = (Aat/ ^a)^n 
with ^=0.66 GeV 2 . The effective confinement size is given by \c a /^N= 1-194 accord- 
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ing to Eq. (|2.1|) and ^c" a =1.74. The confinement size is 5.0 percent larger than the one 
(A(7 /Ajv=1-137) in Ref. 0j. Therefore, other nuclear confinement radii in Ref. H are also 
multiplied by the factor 1.050. To be explicit, they are given for the deuteron, helium, 
carbon, and tin as Xd/^n— 1.066, £d=1.186, A# e /A./v=1.133, £ff e = 1-437, Acf/Ajv=1.159, 
£c=1.557, A5 n /Ajv=l-235, and £s n =2.001, where we assumed that the nuclei are 4 He, 12 C, 
40 Ca, and 118 Sn in the theoretical calculations. All the rescaling parameters are obtained at 
Q 2 =l GeV 2 . Using these parameters, we calculate the rescaling modification. 

In the previous publication [|TT[] , the deuteron structure function F® is assumed to be 



equal to the one for the nucleon in calculating the ratios F^/F®. A special attention should 
be paid for the deuteron if we wish to calculate it in the recombination model. There is 
an explicit nuclear radius factor i? i n Eq. ( |2.2|) and also in the PR evolution equations. 
However, it is well known that the deuteron radius does not follow the rule, R = RqA 1 ^ 3 



with Rq=1.1 fm. We estimate it by the relation R = J 5 < r 2 > /12, where r is the distance 
between the proton and the neutron. Then, the constant Rq becomes i?o=l-56 fm for the 
deuteron. It is fairly large in comparison with the usual nuclear value Rq= 1.1 fm because the 
deuteron is a loosely bound system. It means that the recombination effects in the deuteron 
are much smaller than those in other nuclei. Throughout this paper, we use _Ro=l-56 fm for 
the deuteron. 

With these preparations, we are ready for calculating the F 2 ratios F% / F® ■ First, we 
show the results for the calcium nucleus in Fig. |l|. The dotted curve indicates the distribution 
at Q 2 =l GeV 2 . It should be noted that the convolution integrals of the next-to-leading-order 
(NLO) coefficient functions with the parton distributions are calculated for obtaining the 
i*2- The ratio is evolved to the ones at Q 2 =5 GeV 2 . The dashed curve is obtained by using 
the NLO DGLAP evolution equations, and the solid one is by the PR evolution equations. 
The Q 2 dependence from Q 2 =5 GeV 2 to larger values is rather small as we show it later. 
The experimental data are those of the NMC ||, the Fermilab-E665 |7|], and the SLAC-E139 
pjfl . The theoretical curves are not shown at x < 0.001 because the Q 2 values of the E665 
data are very small and our perturbative calculations cannot be compared with them. We 
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note that there are significant differences between the NMC and E665 data. From Fig. [I], 
we find that our theoretical shadowing is rather small at Q 2 =l GeV 2 in comparison with the 
experimental one. The evolved distribution by the DGLAP shows also smaller shadowing. 
The NMC data can be explained only if the PR evolution equations are used. It is also 
interesting to find that the ratio F^ a /F® >> 1 at large x is explained as a recombination 
effect without the explicit Fermi-motion corrections. We apply the model for other nuclei. 
Our results are compared with experimental data for a small and large nuclei, helium and 
tin, in Figs. ||] and [3] where the EMC data |27[] are also shown. The solid curves indicate our 
theoretical results at Q 2 =5 GeV 2 with the PR evolution. Although the modification tends 
to be small at x ~ 0.7, we find in these figures that our model can explain the x dependence 
reasonably well in the wide x region (0.001 < x < 0.8) for the small, medium, and large size 
nuclei. 

Next, Q 2 dependence of F^ 0, / F® is calculated at the fixed x points, x=0.0085, 0.035, 
and 0.125. Then, the theoretical ratios are compared with the NMC data in Fig. |j. The 
calculated results are shown by the dotted, solid, and dashed curves for :r=0.0085, 0.035, 
and 0.125, respectively. The figure indicates that our model is successful in explaining the 
Q 2 dependence at x=0.035 and 0.12; however, it is not very successful in the small Q 2 region 
at £=0.0085. Namely, the calculated shadowing is too small to explain the experimental one 
in the region Q 2 ~1 GeV 2 . Therefore, the model needs to be studied further whether it is 
in fact impossible to explain the shadowing behavior at small Q 2 . We leave this problem as 
our future research topic. 

We found that our model is successful in explaining the gross properties of the nuclear 
structure functions F^- It is interesting to extend the model to the studies of sea-quark and 
gluon distributions in nuclei. This topic is discussed in the next section. 
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IV. SEA-QUARK AND GLUON DISTRIBUTIONS 



The nuclear modification of sea-quark and gluon distributions is not well investigated. 
Because the F 2 at small x is determined by the sea-quark distributions, the sea-quark shad- 
owing is almost equal to the F 2 shadowing. However, there is no independent experimental 
confirmation of the sea shadowing. There are Fermilab-E772 Drell-Yan data flO)] , which 
measured the ratios, S C /S D , S Ca /S D , and S Fe /S D in the x range 0.04 < x < 0.27. The 
data are shown in Fig. Although the iron data are often quoted in suggesting that there 
is no nuclear modification in the measured x range, it is not very obvious by noting the A 
dependence and the experimental accuracy. The nuclear sea-quark distributions have been 
studied extensively in the context of the pion-excess model ||28|| . However, we do not step 
into the pion-excess problem in this paper. 

We show our theoretical results for the helium, carbon, calcium, and tin nuclei in Fig. 
HI The calculated sea-quark shadowing is not exactly equal to the one for F 2 because the 
gluon distribution also contributes to F 2 through the coefficient function. However, the gluon 
corrections are small so that both shadowing results are almost the same. The ratio becomes 
large at medium x due to the recombination effects. The data for F 2 jF 2 in Figs. [I] and 
indicate antishadowing features at x ~ 0.1. Our model explains this region by the large 
sea-quark ratios S A /S D in Fig. However, we are annoyed that our ratios are too large 
to explain the experimental sea-quark distributions in the region 0.1 < x < 0.2. If there 
were no nuclear modification as suggested by the E772 data, the F 2 antishadowing should 
be interpreted as the valence-quark antishadowing. On the other hand, there seems no 
physics mechanism to produce the valence antishadowing in such a x region [ 2(J . Of course, 
it is possible to produce the antishadowing artificially if the baryon-number conservation is 



imposed with the valence-quark shadowing and the medium-x EMC effect p9 |. However, 
it does not solve the physics problem. We hope that experimentalists measure the sea- 
quark ratios from very small x to medium x accurately for several different nuclei. The 
accurate data should be able to clarify the theoretical problems. At this stage, there is no 
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experimental information on the Q 2 dependence of the sea-quark ratios, of course, except 
for the one implied by the F 2 data at small x. We show the ratios for the calcium nucleus 
at 2=0.0085, 0.035, and 0.125 in Fig. |J The curves in the figure are almost the same as 
those in Fig. [| except at the x=0.125 results where the valence-quark distributions also 
contribute. 

The detailed sea-quark distributions will be investigated experimentally at RHIC (Rel- 
ativistic Heavy Ion Collider). We may wait a few years for the accurate data. In addition, 
flavor distributions of the sea could be an interesting topic. For example, the u/d asymmetry 



in nuclei is important for testing a dynamical aspect of the nuclear parton model |30| . There 
is also a possibility to measure the nuclear u/d asymmetry at Fermilab |JT] . 



Next, the nuclear gluon distributions are calculated. We have already shown the first 



version in Ref. ||16|| : however, the results were not obtained in a consistent manner with the 
F2 structure functions. Here, the model was tested first by comparing various Ff/Ff data. 
We found that it is a fairly successful model in section [nT| , so that it is reasonable to extend 
it to the studies of the sea-quark and gluon distributions. To find the accurate nuclear gluon 
distributions is important not only for establishing the nuclear model at high energies but 
also for other applications to heavy-ion physics. For example, the J/ip suppression may be 
interpreted as a quark-gluon plasma signature; however, the initial gluon distributions in 



heavy nuclei are not known. It would be, at least partially, a local gluon shadowing effect |32 
In order to find the correct interpretation, the information on the precise gluon distributions 
is inevitable. At this stage, most people just assume that nuclear gluon distributions are 
equal to the nucleon's distribution. However, it is known that the recombination mechanism 
produces strong gluon shadowing |I6| , [17| . Despite their importance, it is unfortunate that 
the gluon distributions are not actually measured. There are implicit data on the nuclear 
modification and explicit data taken by the NMC for the ratio G /G . However, the 
data are not accurate enough to enlighten us. 

We show the calculated nuclear modification in Figs. and |8|. Figure |7| shows the x 
dependence for the helium, carbon, calcium, and tin nuclei. The figure indicates that the 
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gluon shadowing is significantly larger than that of quarks, and the ratio becomes large at 
medium x due to the recombination effects. The Q 2 dependence in Fig. [8] shows that the 
gluon shadowing is very large at small Q 2 and it has large Q 2 dependence at small Q 2 . 

It is interesting to test these predictions by future experiments. The RHIC facility is 
appropriate for finding the gluon shadowing ]33J], for example, through the direct-photon 



production process. In order to find accurate sea-quark and gluon distributions, experimen- 
tal information is crucial. 

V. DEUTERON 

The nuclear corrections have been calculated in the deuteron for obtaining the structure- 
function ratios F 2 / F 2 in section |TI. Here, we discuss the correction effects on a sum 



rule. Because there is no fixed target for the neutron, the deuteron or 3 He is usually 
used for measuring neutron structure functions. The nuclear corrections are important in 
extracting the "neutron" structure functions from the deuteron and 3 He data. For example, 
the Gottfried sum rule is given by the difference between F 2 and F 2 [551]: 



Iq= — [Fi(x, Q 2 ) - F?(x, Q 2 )} = - + - / dx [u(x, Q 2 ) - d(x, Q 2 )] . (5.1) 
Jo x 3 3 Jo 

The NMC used the measured proton and deuteron structure functions, F 2 and F 2 D , for 



evaluating the difference [36 



F 2 p -F 2 n = 2F 2 D l i + fJ / F J i , F?/FZ = 2F 2 D /F!-1. (5.2) 
These equations, of course, assume no nuclear correction. As it was discussed in section 



III], the nuclear structure functions are modified in the whole x region. In particular, it is 
known that the small-x part, namely the deuteron shadowing, is important for obtaining 
the "correct" sum in Eq. ( |5.1| ). 

The F 2 / F 2 ratios are measured by various groups. We compare our theoretical results 



with the data by the NMC and the Fermilab-E665 PS] in Fig. PI. The dotted curve 
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is the MRS-R2 distribution (u ^ d) for the F?/F% ratio at Q 2 =l GeV 2 , and the dashed 
curve is the flavor symmetric one which is given by the modified MRS-R2 so as to have 
u = d = s = Sea/6 at Q 2 =l GeV 2 . The calculated distribution 2F 2 D /F 2 D - 1 at Q 2 =5 
GeV 2 is shown by the solid curve, which includes the nuclear corrections. Because the E665 
small-x data are taken at very small Q 2 where the perturbative QCD would not work, the 
theoretical curves are not plotted at small x (<0.001). 

Because the MRS-R2 was created so as to explain the NMC data as well as many other 
data, the dotted curve is in agreement with the NMC ratios. The flavor symmetric distri- 
bution, which is shown by the dashed curve, is well below the NMC data at small x where 
sea-quark contributions are significant. The solid curve includes the nuclear modification 
due to the rescaling and recombinations, and it also contains Q 2 evolution effects from 
Q 2 =l GeV 2 to 5 GeV 2 . As a result, the solid curve is below the NMC data in the x region 
(0.0015 < x < 0.675). It is mainly because our model used the flavor symmetric input in 
order to see the nuclear modification effects. However, the difference is partly attributed to 
the incomplete analyses by the NMC and MRS-R2 in the sense that the nuclear corrections 
are completely neglected. In order to illustrate the significance of the corrections, the nuclear 
contribution to the Gottfried sum is estimated by 

Ia = lS+f 1 _ ^2[F 2 D (x)-F 2 N (x)]. (5.3) 

The Iq is the sum without the nuclear correction, and the second term (= 5Ig) is the 
correction. The integral is evaluated for the structure functions at Q 2 =5 GeV 2 , and we 
obtain 5Ig = —0.013 for the integral from x m i n = 0.004. It is the same order of magnitude 



as the other recent shadowing estimates P9fl . This kind of correction should be taking into 
account in discussing the Gottfried sum rule and also other sum rules with neutron structure 
functions. 
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VI. SUMMARY 



We have investigated the parton model with the Q 2 rescaling and parton-recombination 
mechanisms for describing the parton distributions in nuclei. The nuclear modification of the 
x distributions is calculated at Q 2 =l GeV 2 . If they are evolved to larger Q 2 , for example 
Q 2 =5 GeV 2 , by the parton-recombination evolution equations, they become comparable 
with the experimental F 2 shadowing. 

The model was then extended to the studies of sea-quark and gluon distributions in 
nuclei. We showed the x and Q 2 dependent results of the ratios Sa/S d and Ga/G d . They 
indicate large nuclear shadowing at small x and large ratios in the medium-x region. These 
predictions could be tested by the experiments at RHIC. The understanding of nuclear 
parton distributions is very important for heavy-ion physics studies, for example, in finding 
a signature of the quark-gluon plasma. 

We also indicated in our model that the nuclear effects should be taking into account 
properly in order to discuss the sum rules with neutron structure functions, which are 
extracted from the deuteron data. 
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FIG. 1. Theoretical results for the ratio F^/F® are compared with the NMC, Fermilab-E665, 
and SLAC-E139 data. The dotted curve is the initial ratio at Q 2 =l GeV 2 . The distributions are 
evolved to Q 2 =5 GeV 2 , and the dashed and solid curves are obtained by the DGLAP and PR 
evolution equations, respectively. 
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FIG. 2. Nuclear modification of Fi is shown for the helium nucleus. The theoretical result is 
shown by the solid curve at Q 2 =5 GeV 2 . 
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FIG. 3. The modification is calculated for the tin nucleus. The theoretical result is shown by 
the solid curve at Q 2 =5 GeV 2 . The E665 data are also shown even though they are taken for a 
different nucleus, xenon. 
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FIG. 4. Calculated Q 2 dependence of the ratio F^jF® is compared with the NMC data. The 
dotted, solid, and dashed curves are obtained at x=0.0085, 0.035, and 0.125, respectively. 
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FIG. 5. Calculated x dependence of the sea-quark modification is shown by the dotted, dashed, 
solid, and dot-dashed curves for the ratios S He /S D , S C /S D , S Ca /S D , and S Sn /S D at Q 2 =5 GeV 2 . 
They are compared with the Fermilab-E772 data for the carbon, calcium, and iron. 
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FIG. 6. Calculated Q 2 dependence of the sea-quark modification is shown for the cal- 
cium-deuteron ratio S Ca /S D . The dotted, solid, and dashed curves are obtained at x=0.0085, 
0.035, and 0.125, respectively. 
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FIG. 8. Calculated Q 2 dependence of gluon modification is shown for the calcium-deuteron 
ratio G Ca /G D . The dotted, solid, and dashed curves are obtained at £=0.0085, 0.035, and 0.125. 
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FIG. 9. The ratio of the neutron and proton structure functions is shown. The dotted curve 
is the MRS-R2 at Q 2 =l GeV 2 . The flavor symmetric one is shown by the dashed curve at Q 2 =l 
GeV 2 . The solid curve is the evolution result at Q 2 =5 GeV 2 , and it is calculated by 2F® /F% — 1 
with the nuclear corrections in the deuteron. 
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